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We theoretically investigate phosphorene zigzag nanoribbons as a platform for constriction engi¬ 
neering. In the presence of a constriction at one of the edges, quantnm confinement of edge protected 
states reveals conductance peaks, if the edge is uncoupled to the other. If the constriction is narrow 
enough to promote coupling between edges, it gives rise to Fano-like as well as anti-resonances in 
the transmission spectrum. These effects are shown to mimic an atomic chain like behavior in a two 
dimensional atomic crystal. 
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I. INTRODUCTION 

Low dimensional systems have attracted the attention 
at least in the past fifty years, since the development of 
semiconductor epitaxial growth and deposition of metal¬ 
lic thin films^. The early scenario, back in the 1960’s, 
as promising as it appeared, has evolved into a main¬ 
stream interest in condensed matter physics due to land¬ 
mark discoveries in the late 1970s and early 1980s, like 
the quantum Hall effect^ and the conductive polymers^, 
respectively 2D and ID systems. The subsequent discov¬ 
ery of new carbon allotropes, showing stable structures 
either in OD (fullerenes), ID (carbon nanotubes) and 2D 
(graphene) consolidated this scenario in an exciting re¬ 
search field^. The isolation of strictly one atom thick lay¬ 
ers in the first years of the present century opened a wider 
window for both basic physics and device applications®. 
These new disruptive research efforts, initially impulsed 
by graphene, are nowadays detaching from carbon based 
roadmaps, as also envisaged, at the beginning of the 
graphene boost, by Novoselov, Geim and coworkers®. 

A very recent 2D atomic crystal of black 
phosphorousnamely phosphorene, is a promis¬ 
ing system in which 2D properties together with strictly 
ID chain behavior are present in different energy win¬ 
dows, hence allowing to a same device to be tuned from 
a ID to a 2D system by simply tuning the Fermi energy. 
In the present work we focus on the ID energy window, 
created by effective doubly degenerate band (in the rele¬ 
vant energy scale) associated to states strongly localized 
at the zigzag edges^^’^^ of phosphorene nanoribbons, 
whose properties are explored using a new strictly one 
dimensional resonant tunnelling device. 

The double barrier resonant tunnelling device^’^®’^^, 
conceived here as an atomically precise segmentation at 
one of the edges, shows unusual geometry, since the di¬ 
rection of the barriers is perpendicular to the well and 
contact regions^®’^®. Among our findings we show that 
for a thin barrier case (constriction with narrow step 
from the upper zigzag edge), the resonant tunnelling 
permits a spectroscopy of the band structure of phos¬ 


phorene nanoribbons in this energy window. Further¬ 
more, progressive widening of the barriers (enhancing the 
step width of the constriction) therefore nearing the con¬ 
striction to the other edge leads to edge coupling effects 
featuring resonances with Fano line shapes^^^^® reveal¬ 
ing also a new discrete/continuum states coupling sys¬ 
tem. For this latter coupled edge system, the transmis¬ 
sion probability characteristics turn out to present clear 
features of both (i) the actual finite confining segment 
coupled to an infinite (not segmented) edge and, (ii) the 
properties of an infinite narrow nanoribbon with strongly 
coupled edges. These results are independent of the area 
of the device region, solely from the segmented region 
length and distance between the edges, revealing an ef¬ 
fective chain-like behavior of the nanoribbon’s edges. 

In what follows, we initially discuss the “bulk” elec¬ 
tronic properties of a phosphorene nanoribbon, present¬ 
ing the model calculation framework, as well as the ef¬ 
fects of edge coupling on the conductance of these infinite 
zigzag ribbons that are essential to understand the reso¬ 
nant tunnelling spectra. Next, the geometry of the actual 
investigated segmented device is presented, introducing 
the resonant tunnelling effects. In the sequence, the core 
of the results is devoted to explore the segmented edge 
device resonant tunnelling behavior, showing how the de¬ 
fects on the edge may actually enrich the scenario instead 
of solely washing out the announced effects, giving fur¬ 
ther evidences that the phenomenon is restricted to the 
atoms at the very edge. Finally, the conclusions suggest a 
bridge between the present ID systems embedded in a 2D 
crystal and ongoing research on isolated atomic chains. 


II. PHOSPHORENE ZIGZAG NANORRIBONS 
AND MODEL CALCULATION: EDGE 
COUPLING EFFECTS 

The essential atomistic aspects of the structures inves¬ 
tigated are depicted in Fig. 1. Fig. 1(a) shows a segment 
of an infinite zigzag edged phosphorene nanoribbon of 
width Nz = 8, the number of zigzag chains along the rib- 
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bon in this case. The tight binding hopping parameters 
considered, as discussed below, are indicated in Fig. 1(b). 

The quite complex electronic structure of phosphorene, 
already at energy ranges rather close to the Fermi en¬ 
ergy, hinders a wider use of single orbital tight-binding 
models in chasing the alluded electronic and transport 
properties of systems based on this new material. Never¬ 
theless, the use of such model is well validated, by means 
of comparisons with first principle electronic structure 
calculations^*^, for the very energy window of interest 
around the gap, namely a double central band. This cen¬ 
tral band for zigzag phosphorene nanorribons has been 
predicted for phosphorene*^^’^^ and is absent in Graphene. 
We use here the same tight-binding parametrization 
for phosphorene proposed by Rudenko^*^ considering a 
Hamiltonian H = where Ci (c|) is the cre¬ 

ation (annihilation) electronic operator at site i and tij 
is the hopping integral between sites i and j. In this 
model, to characterize the low energy electronic proper¬ 
ties five hopping integrals are required^*^: ti = —1.220 
eV, t 2 = 3.665 eV, tz = -0.205 eV, ^4 =-0.105 eV and 
= -0.055 eV . The transmission T = Tr[riG’TijG“] is 
calculated using the recursive Green’s function^^ O'" = 
[E + ir] — H — T,]^ — in the phosphorene lat¬ 

tice representation. The left and right contacts broad¬ 
ening function and the self¬ 

energy of contact ^l{r) are calculated recursively for the 
semi-infinite zigzag phosphorene nanoribbons^^; other 
electronic properties such as the local density of states 
(LDOS) Pa = —Im[G’'(ri, 7y, i?)]/7r are also calculated. 

The electronic and transport properties of a host 
zizgzag nanoribbon in which a finite segment will be lat¬ 
ter tailored in are also summarized in Fig. 1. The inset 
in Fig. 1(c) depicts the energy window of interest, show¬ 
ing the top(bottom) of the valence(conduction) band and 
an effectively degenerate central band*^^. These central 
bands present cosine like dispersions characteristic for 
ID systems^^. Indeed, the degeneracy comes from the 
fact that the width of the ribbon here is Nz = 60, which 
guaranties that the two edges are effectively uncoupled^*^. 
Hence this width will be chosen for the host ribbon where 
the constriction will be introduced. 

The effect of edges coupling on the band structure can 
be followed in the main part of Fig. 1(c), showing a zoom 
of the central band energy range. Having in mind the un¬ 
coupled limit of Nz = 60 (red curve), lifting of the central 
band effective degeneracy starts (in the present relevant 
energy scale) for Nz = 15 (black curve) at the center 
of the Brillouin zone with an approximately symmetric 
splitting of slightly deformed cosine-like bands. Indeed 
an incipient overlap of wave function in this situation is 
illustrated in Fig. 1(e), with a noticeable amplitude of 
the wave function in the atomic sites well inside the rib¬ 
bon (It should be noticed that Nz corresponds to the 
number of zigzag chains, hence for Nz =15, there will 
be 30 atomic sites in the unit cell). For extremely thin 
ribbons, Nz = 7, the splitting becomes of the order of 
the uncoupled band widths (blue curve). More striking 



FIG. 1. (Color online) (a) Illustration of a zigzag phos¬ 
phorene nanoribbon of width N^, where Nz is the number 
zigzag chains, (b) Hopping parameters used in the four band 
model^**. (c) Structure of the central band of edge states 
of nanoribbons with different widths. The inset includes the 
band structure of the bottom and top of the conduction and 
valence bands for the Nz = 60 case, (d) Transmission prob¬ 
abilities for the edge states depicted in (c). (e) Probability 
amplitude of edge states of the Nz = 15 case for different ka 
values. The inset reveals the stronger non-zero overlapping 
of states of the two edges at the bottom of the central band, 
ka = 0. 


is the drastic change in the shape of one of the bands, 
showing a local maximum at the center at Ka = 0 and 
a minimum at Ka = 1. 

The consequences of the edge coupling on the transmis¬ 
sion probabilities along the edges are qualitatively signifi¬ 
cant as can be seen coming back to Fig. 1(d). Degenerate 
bands sum up to a plateau of T = 2 (red curve). A slight 
lifting of the degeneracy breaks the lower threshold of the 
plateau introducing a T = 1 step, an energy range where 
there is only one conduction channeh*^. However, ex¬ 
treme coupling leads to a T = 3 plateau for — 1 < fca < 1, 
where the different states corresponding to the same en¬ 
ergy in this band are added to the channel associated to 
the other cosine-like band. 
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FIG. 2. (Color online) (a) Schematic picture of a constric¬ 
tion characterized by the parameters L =10, actually used 
throughout the work, and ruz = 4, here only for the sake of 
illustration, (b) Transmission probabilities, as a function of 
the Fermi energy, at the energy range corresponding to the 
central band, including the bottom and top of the conduction 
and valence bands. Three different situations are depicted: 
two constriction defined at nanoribbons of width Nz =60, 
both with I/=10, and rriz = 57 (red) and rriz = 15 (green) 
and a zigzag nanoribbon without any constriction, Nz = ruz 
= 57 (black line). 


III. SEGMENTED NANORRIBONS: 

RESONANT TUNNELLING IN ID EFFECTIVE 
CHAINS STRUCTURES 

In the energy energy range of the edge states band, 
the “bulk” of the nanoribbon acts mainly as “in plane” 
substrate for the two dimensional channels at the edges. 
This condition, evidenced by the electronic band struc¬ 
ture discussed in the previous section raises the ques¬ 
tion of a mean to observe experimentally those effective 
one-dimensional chains embedded in the rather complex 
phosphorene crystalline structure. In order to test our 
hypothesis we propose the segmented nanoribbon struc¬ 
ture (constriction) illustrated in Fig. 2(a). The segment 
of a thinner region of the nanoribbon of width mz, also 
given in number of zigzag chains, is defined by a length 
L in units of atoms removed along one zigzag direction. 
One essential parameter is the step width between the 
semi infinite upper edges and the central segment, which 
is simply defined as Nz — mz and, as will be seen below, 
defines the barrier thickness in the resonant tunnelling. 

shows the transmission probabilities through two con¬ 
strictions of lenght L = 10 with Nz — mz = 3 and 
Nz — mz = 45 step widths, compared to a bare N^ = 57 
nanorribon, as a function of energy. In order to avoid 
any coupling between the edge states the width of the 


nanoribbon in the contacts is also N^ = 60. Transmission 
plateaus above (below) the edge states band are shown, 
for the sake of completeness, since these structures are 
of entirely different character than the resonances in the 
central band we will be focusing on. These transmis¬ 
sion plateaus are due to the lateral confinement in a 
nanoribbon. This is confirmed by the green curve for 
mz = 15, a deep step leading to a large shifting of the 
valence and conduction bands transmission plateaus, ev¬ 
idencing also the well known Fabry-Perot oscillations^^ 
due to the geometrical modulation; these effects are al¬ 
ready well known for graphene and square lattice nanor- 
ribons with constrictions^"^^^®. 

On the other hand, the edge states, observed here in 
the energy range from 0.3 eV to 0, drops to T[E) =1 with 
resonant peaks on top for wide constrictions and anti- 
ressonances for narrow constrictions. Fig. 3(a) presents 
a closer look at this energy range; for thin armchair 
steps {Nz — mz = 3) the red curve shows a group of 10 
peaks, these resonances get thinner as the step increases 
{Nz — mz = 5). Further diminishing ruz, the barriers 
to the upper edge contacts become too large, but now 
the coupling to the lower edge becomes relevant. For 
mz <15 transmission shows asymmetric Fano-like reso¬ 
nances at the low energy side and sharp anti-resonances 
at higher energies within the central band as will be dis¬ 
cussed below. It should be recalled here that only for 
extremely thin nanoribbons the strong coupling between 
the edges widens the central band, however, it is not ob¬ 
served in Fig. 3(b) a plateau enlargement for = 8 
because the edges of the left and right contact are not 
coupled (keeping the central band of the contacts unal¬ 
tered). 


In our constriction the role of both channels (discrete 
states and continuum at the upper and lower edges re¬ 
spectively) can be made explicit by picturing the local 
density of states (LDOS) in Fig. 4, at the energy values 
indicated by arrows 1, 2, 3 and 4 in Fig. 3. 

In Fig. 4(a), looking at the LDOS corresponding to a 
transmission peak pointed out by the arrow 1, it is clear 
that the higher values of LDOS appear on the edges, 
the confined state at the constriction in the upper edge 
clearly stands out. It should be noticed that this LDOS 
is slightly asymmetric, since the structure with L equal 
to an even number of atoms is asymmetric (see Fig. 2(a)). 
This asymmetry leads to a resonance peak T < 1^® su¬ 
perposed to the background plateau. L equals to an odd 
number of atoms restore complete symmetry, leading to 
higher resonances, T « 1 (not shown here). A less in¬ 
tense LDOS along the lower edge corresponding to the 
T = 1 plateau can also be observed. The LDOS along the 
upper edge outside the constriction region is less intense 
in the figure scale, due to the prominence of the confined 
state. Far from a resonance, actually between two reso¬ 
nances, situation pointed out by arrow 2 in Fig. 3(a), the 
LDOS in the confined part of the upper edge is strongly 
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FIG. 3. (Color online) (a) Transmission probabilities, as a 
fnnction of the Fermi energy, at the energy range of the central 
band corresponding to nanoribbons with constrictions: red 
and green correspond to the cases in Fig. 2(b): =60, all 

with L =10 and = 57 ( shallow constriction) and = 15 
(deep constriction), respectively. The black cnrve corresponds 
to a constriction with = 55 (intermediate depth), (b) 
Transmission probabilities for very deep constrictions, = 
8 and different lengths: Z/=10 (upper panel) and L =30 (lower 
panel). 


suppressed, enhancing the contribution along the entire 
lower edge and the upper edge contact sections (i.e., out¬ 
side the confining region), Fig. 4(b). 

A quite different situation is depicted in Fig. 4(c)-(d), 
where the confined sate in the constriction is decoupled 
from the upper edge with a varible coupling to the lower 
edge. In Fig. 4(c), corresponding to the anti-resonance 
labelled as 3 in Fig. 3, we observe a faint coupling to the 
lower edge. 


The LDOS plots reveal the unique character of the res¬ 
onances observed in Fig. 3: resonant tunnelling through 
confined edge states in the constriction defined by an 
armchair-like step double barrier structure. It should be 
noticed that the number of resonances is identical to the 
number of atoms removed along the segment that define 
the length of the constriction, L — 10, indicating that 
the transmission shows a spectroscopy of the ID states 
at the edge of the constriction. Indeed, increasing the 



FIG. 4. (Golor online) LDOS at the energies pointed out by 
arrows, labeled 1, 2, 3 and 4, in Fig. 3 for constrictions L=10 
long. The upper panel is for shallow constrictions, = 60 
and rriz = 57: (a) at a resonant energy, corresponding to 
arrow 1. (b) Off resonant energy, corresponding to arrow 2. 
The lower panels are for deep constrictions: (c) at the anti 
resonance indicated by arrow 3 (m^ = 15); (d) at the anti 
resonance highlighted by arrow 4 (m^ =8). 


length of the constriction will increase at the same pro¬ 
portion the number of the resonances (not shown here). 
The fact that the resonances are insensible to the con¬ 
striction width (hence the area for a fixed length) is a 
further indication that we are dealing with a strictly ID 
effect at the edges: albeit the underlying 2D crystal, the 
behavior revealed here is the one of an effective atomic 
chain. Here we should notice that the only signature 
of the underlying 2D crystal is given by the resonances 
widths. Recalling Fig. 1(e), the resonances near the bot¬ 
tom of the central band, Ka « 0, correspond to states 
which deeper penetrate in the bulk, hence the barriers 
defined by the device steps are less effective than for res¬ 
onances at higher energies. Two dimensional structured 
systems, like nanorribons, present transmission proba¬ 
bilities with multi channel contributions that are mixed 
by the geometrical changes along the structure^®’^^. On 
the other hand, one dimensional systems present sin¬ 
gle channel transmission probabilities, that are described 
by s-like orbital chain models. In the present case, in¬ 
deed, the positions of the transmission resonances are 
shown in Fig. 5 to be well reproduced by a simple one 
dimensional double barrier quantum well modelled by 
a chain of s-like orbitals^^’^^. The nearest neighbor 
hopping (til)) of the one dimensional chain, shown in 
Fig. 5(a), is calculated by |ti_D| = Aif/4 = 0.0775 eV, 
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FIG. 5. (Color online)Comparison between the transmission 
probabilities at the central band energy ranges for shallow 
constrictions, niz = 57, and equivalent s-like orbitals chain 
toy model, (a) Representation of the one ID double barrier 
quantum well structure at the upper edge, the central blue 
sites represent the quantum well at the constriction, while the 
red ones are for the barriers (armchair steps). The left and 
right contacts at the upper edges are also represented by blue 
sites., (b) Transmission probabilities for a L =10 constric¬ 
tion: 4 band tight binding model (red) and the toy model 
(blue), (c) The same as (b) but for a longer constriction, L 
=15. 


where AE is the edge states band width obtained from 
the red curve in Fig. 1(c). The transmission peaks of 
the one dimensional quantum well with L = 10 and 
L = 15 atomic sites clearly reproduce the position of 
the resonaces observed for phosphorene constrictions of 
the same length, see Fig. 5(b)-(c). The resonant peaks 
appear at the energies of the infinite square well energies 
En = ‘2tiD + 2ti£)COs(n7r/(L -|- 1)) for n = 1, 2,..., L 
To recap, the edge states at opposite edge provide one 
dimensional electronic transport channels embedded in a 
two dimensional material and the conductance observed 
can be understood from simple model. The resonant 
peaks on top of the T = 1 plateau resemble the con¬ 
ductance of two parallel and independent channels, as 
shown in Fig. 6(a). The lower edge provides a contin¬ 
uous channel of T = 1 while the upper edge presents a 
quantum well with tunnelling coefficients 7^(7?) accross 
the left(right) barrier defined by the vertical edges in the 
figure. Hence, resonant tunneling becomes only possi¬ 
ble when the energy matches the energy of the bound 
states in the well. When this situation is not fulfilled 
the upper channel is closed and the transmission of the 
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FIG. 6. (Color online) Schematic representations of the main 
transmission probability conditions. The grey areas represent 
the nanoribbons with constrictions, with different edge states 
behaviors highlighted by red and green lines, (a) Resonant 
transmission at the upper edge summed up to the continuous 
transmission at the lower edge, hence T =2. The resonant 
coupling of the confined state at the constriction with the 
left(right) contact is represented by 7i(7ij). (a) Off resonance 
negligible transmission at the upper edge with the continuous 
transmission at the lower edge, T =1. (c) Anti-resonance in 
the transmission due to the strong coupling, represented by 
Fs, between the confined state at the deep constriction and 
the lower edge, now depicted as a dashed line indicating the 
absence of transmission, T = 0. 

whole system is T = 1, see Fig. 6(b). A deeper step on 
the constriction, on the other hand, would lead to quasi 
bound states, i.e., wider barriers, connected to the up¬ 
per edge contacts, but with a significant coupling, 73, 
to the bottom edge. Fig. 6(c), leading to Fano-like reso¬ 
nances and antiresonances in the transmission^^’^"^. This 
situation corresponds to quasi -bound state coupled to 
a continuum, leadind to fano-like asymmetric resonances 
and antiresonances in the transmission. 

The three different line-shapes can be described by a 
single expression^®: 

T{e) =1 t, p ^ (1) 

where td is the direct transmission without the pres¬ 
ence of a scattering region, e = {E — Eji)/T, {Er is the 
energy of the resonant discrete state and F its line width) 
and q is the Fano asymmetry factor,which represents the 
ratio of the resonant tunneling channel to the channel 
due to the continuum (here represented by the lower edge 
channel). When g —>■ 00 and —>■ 0 (no continuum chan¬ 
nel available) a a resonance peak develops; for g « ±1 
(both channels are relevant), an asymmetric line shape 
is revealed, while for g = 0 an antiresonance appears. 
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Recalling the framework of the present work, edge con¬ 
fined states are supported only by zigzag edges, being 
absent in armchair or bearded edges^^d2 Therefore, in¬ 
troducing perturbations to a zigzag edge, like edge va¬ 
cancies, would locally destroy these ID states. The con¬ 
sequences of these perturbations are very relevant in the 
present early stage of phosphorene experimental devel¬ 
opment, in which only preliminary steps toward design 
and realization of effective devices out of the bulk in the 
nanoscale have been reported so far.^®. Hence, effects of 
the presence of disorder at the edges have to be consid¬ 
ered. 

In Fig. 7 we present the transmission probability as a 
function of energy as well as LDOS associated to selected 
resonances in the presence of vacancies. 

Defects are normally seen as mechanisms that hin¬ 
der the observation of transport properties associated 
to shape modulation of nanoscopic low dimensional sys¬ 
tems. Indeed, the resonance spectra are also dramati¬ 
cally modified in the present case. However, the issue 
can be seen from an entirely different point of view. The 
vacancies change locally the character of the edge, in¬ 
troducing actually small barriers, further dividing the 
system into smaller segments. The system chosen here 
is a device with a vacancy located at the upper edge of 
the left contact near the central segment (quantum well), 
the exact position is marked by the arrow in Fig. 7(d). 
What can be observed in the transmission probabilities 
in Fig. 7(a)-(c) is that the resonances of the well at the 
contact defined to the left by a barrier due to the va¬ 
cancy, couple to some of the states of the original well 
given by the central segment. Those couplings are iden¬ 
tified by the split peaks clearly seen in Fig. 7(b). If the 
barriers defining the central segment (well) are widened. 
Fig. 7(a), the splitting diminishes. Fig. 7(c) depicts the 
same device in absence of the vacancy as a guide for iden¬ 
tifying the couplings. The character of the ID confined 
states in the presence of a vacancy is illustrated in the 
LDOS, Fig. 7(d)-(e), for the resonances corresponding to 
arrows I and 2 in Fig. 7(b), respectively. State 1 corre¬ 
sponds to a double-well (vacancy-upper edge contact-left 
step-central segment-right step) along the upper edge, 
exhibiting LDOS at both wells, while state 2 is confined 
mainly to the central quantum well. 

Having in mind the previous discussion, while a res¬ 
onant tunnelling spectroscopy would become rather in¬ 
volved with the presence of defect induced barriers, scan¬ 
ning probe microscopy continues a way to reveal the edge 
quantum confinement. 


IV. CONCLUSIONS 

It is inevitable to compare our results with those ob¬ 
tained for graphene constrictions: both zigzag graphene 
and phosphorene nanoribbons support edge states, how- 
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FIG. 7. (Color online) Edge states transmission probabilities 
for a nanoribbon Nz, = 60 wide, with a constriction L =10 
long, adding an edge vacancy at the left contact: (a) rriz = 
55 and (b) tHz =57. (c) The system without the vacancy 
for the sake of comparison, highlighting the splitting of the 
resonances due to the presence of a vacancy, (d) LDOS at 
energy indicated by arrow 1 in (b). (e) LDOS at the energy 
indicated by arrow 2 in (b). 


ever, their signatures on the electronic transport proper¬ 
ties are completely different. First, edge states in zigzag 
graphene nanoribbons are sublattice polarized, so one 
single edge do not contribute to the electron transport 
properties. The graphene edge states channel is orig¬ 
inated by the overlapping of edge states on opposite 
edges^®, contrary to what we observe here where a sin¬ 
gle phosphorene edge provides an independent transport 
channel. Second, localized states in graphene junctions 
manifest as antiresonances of zero conductance, these 
states localize over the junction^®^^^, while the localized 
states of phosphorene constrictions strictly on grooved 
zigzag edge and appears as peaks, asymmetric Fano line 
shapes or dips in the conductance. In summary, we pro¬ 
pose phosphorene zigzag nanorribons as a platform for 
constriction (segment) engineering. In the presence of 
an engraved segment at the upper edge, quantum con¬ 
finement of edge protected states reveals resonant tun¬ 
nelling transmission peaks if the upper edge of the host 
nanoribbon is uncoupled to the lower edge. Coupling 
between edges in thin constrictions give rise to Fano- 
like and antiresonances in the transmission spectrum of 
the system. One could envisage to look these effects 
by means of transport measurements as well as scan- 
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ning probe microscopy^^. The energy scale given by the 
resonance spacing is of the order of 10 meV for constric¬ 
tion lengths of L =30 (not shown here), corresponding to 
^ 5.05 nm and contacts ^ 13.3 nm wide, a benchmark for 
experimental efforts, recalling that defects may lead to 
more complex spectra without washing out the main fea¬ 
tures. Other resonant tunnelling mechanisms have been 
observed in phosphorene nanoribbons with vacancies^®, 
defects®® and transverse electric helds^®. It is important 
to reinforce that these mechanisms involve states at the 
interior of the nanoribbon, whereas the effect shown here 
requires one dimensional states localized at the edges. 

Concomitant to the development of the fascinating 
physics of 2D materials, new extreme ID systems, namely 
isolated atomic chains, either based on Carbon^® or 
metallic elements^^, have been obtained and character¬ 


ized, with their properties and possible applications the¬ 
oretically investigated. The present results suggest a way 
were effective atomic chains come out from the edges of 
a 2D crystal. 
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